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SUS? 1 ARY 


This study was carried out in an effort to cetermine the 
longitudinal stability derivatives of a Ryan Navion airplane using dynamic 
flight testing methods. This was cone by recording the transient response 
of the airplane to elevator forcing functions in the form of steps or 
pulses. Data reduction was carried out using the equations of motion 
method and the derivative method. Due to partial failure of the instru- 
mentation the results from the equations of motion method were disappointing. 
More successful results were obtained using the derivative method in con- 
junction with the useful data acquired from the equation of motion method. 

Concurrent with this report an investigation was carried out to 
find the stability derivatives of the same airplane from both theoretical 
considerations anc from steady state flight testing. This investigation 
is described in Reference 7. An effort has been made to correlate the 
flight test results obtained from the dynamic and steady state methods and 
predictions made from theory. This correlation appears in Appendix B. 

Reasonable values of CL were obtained from dynamic testing 
which were consistent with speedc,and position of center of gravity changes. 

Dynamic testing also gave reasonable values for the damping 
derivatives which were consistent at the same speed, but which indicatec 
an increase in the damping derivatives with increase in speed. 

The proper changes in values for the stick fixed maneuver margin 
were found as the position of the center of gravity was changed, however 
a decrease in maneuver margin was found as speed was increased. 

The elevator power was found with good consistency as the center 
of gravity was changed, but decreased in value as the speed was increased. 


a4 





All values of Cm s were somewhat high as comoared with theoretical values. 
An attempt was mace to explain the changes in Cms and stick 
fixed maneuver margin that occurred with speed changes in terms of aircraft 


non-linearities. 
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LIST OF SYMBOLS 


Oy 


(S s 2 Cm partial derivative of the moment coefficient with respect 
< 2% to angle of attack. 


steady state Lift coefficient. 


g 








a 5 re partial derivative of the moment coefficient with respect 
Md ddx todana=vnc, 


Gg. = 2 Cm partial derivative of the moment coefficient with respect 
tee Me ere 


Cot: 2Ge elevator power = partial derivative of the moment 
a0 coefficient with respect to elevator angle. 


On ao slope of the lift curve 
x. 
ee (Pitching moment, ft. lbs.) 
Be 2 {ae oO 
Coes perturbation angle of attack, radians. 
C= perturbation angle of pitch (positive nose up) radians. 


6 = IEG ty, -6 radians. 


dO = vo radians, 


=e, 


7 = perturbation normal acceleration \positive for push over 
maneuvers ). 
Jn 


/ 
nde = [nde 
S = perturbation elevator angle, radians. 


[$= "d%e-#SSdt 


fF. = stick force, lbs. 
Jf. = rate of change of downwash with angle of attack. 


VV = velocity, feet per second. 


V 
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SS (Dh Om cn 


oe 
a \\ 


Dihedral 


mass of the airplane, slugs. 


“ps e 


Me sy seconds. 


moment of inertia of the airplane about the Y-axis, slug me 


2 ky a2 


Lh » radius of gyration, feet. 


mean aerodynamic chore, feet. 
wing area, square feet. 


tail length, distance from the center of gravity to the aero- 
dynamic center of the tail, feet. 


134.2 rt.” = wine area. 
43 rt” = horizontal tail area. 
15.04 rt’ = elevator area. 


6.94 = aspect ratio of wing 


= 3.98 = aspect ratio of horizontal tail. 


33.33 ft = wing span 


13.17 ft = tail span 
9-7 ft = mean aerodynamic chord of wing. 
3.54 ft = mean aerodynamic chord of horizontal tail. 


15.94 ft = tail length, distance between center of gravity and 
aerodynamic center of the tail. 


‘= «24 = position of wing aerodynamic center in % 


- 2° at root = incidence of wing 


lav tip 


-3° = incidence of horizontal stabilizer. 


0 
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INTRODUCTION 

Dynamic testing of aircraft for stability parameters has received 
increasec emphasis with the advent of guided missiles, sonic, and super- 
sonic aircraft. Some of the major efforts in this field have been dynanic 
testing of a Be25 by Marshall E. Mullins at tdwards Air Force Base, (Ref. 
4)3 the dynamic testing of an F-30A by R. C. Kidder at Cornell Aeronautical 
Laboratories, Inc. (Ref. 5)3 and guided missile testing by The Naval Air 
Missile Test Center, Pt. Mugu, California. In all previous cases no effort 
has been made to compare the results with those obtained by steady state 
flight test methods for the same aircraft. This report is part of a 
triple project on a Ryan Navion in which dynamic results are to be conpared 
with theoretical and steady state results. The actual comparisons are 
included in Apvendix B. 

The dynamic testing with which this report 1s concerned consists 
of analysis of the transient response of the aircraft to an elevator forcing 
function. Using various forms of the equations of motion the stability 
parameters for the airplane were reduced from the data by the method of 


averages, least squares, and by the use of graphical plots. 
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THEORY AND ANALYSIS 


The normal procedure in the determination of aircraft transient 
motions is to solve the differential equations which cescribe the aircraft 
motion. In doing this the constant coefficients in the equations of motion 
are estimated from theoretical considerations or measured by steady state 
test methods. 

The problem dealt with in this report is the inverse of that 
stated above; ie., the determination of the coefficients in the equations 
of motion if the time history of the aircraft transient motions are known. 
The response history of the aircraft can be displayed as a function of 
frequency as well as time. In theory the response in either the time or 
frequency domain can be analyzed to determine the constant coefficients 
in the equations of motion. This report deals with the determination of 
the constant coefficients of the longitudinal equations of motion from 


transient responses cisplayed in the time domain. 


Development of Equations: 
Throughout this report it is assumed that the incremental changes 
in forward velocity are negligible. With this assumption the equations 


of motion as developed in Ref. 1 are: 


“ia + da -d® =o (1) 


Cm X + Cmg yd & a Cm 92? - hae Cg =O (2) 
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If elevator lift is consicered important, its effect can be included in 
the lift equation since: 

iter. 

5 4, (3) 





Substitution in equation (1) yields: 





“st + dod - dO - — 5. =O 


oe (4) 


Solving equation (4) for dol and substitution in equation (2) yields: 


(Cn CrragC eg) of + (Cry Crnyg)d@ -hd® + Cong (14 cops) $20 


(5) 


Normal acceleration, Nn, anc elevator angle, 5b, may be substituted for « 
since: 
Cc Cms 


oes nN + (6) 
CL Cie 








Substitution of this expression for & into equation (5) yields; 


Gere Crt Cm Lae 
—- Sa} n + aa Cm, ,) 6 + Con, (i*%,)é hd@=O (7) 
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This equation if integrated once with respect to Ye. yields: 








i Cm, Ga u 
Cee SG) Jd + (Coat Gr) 8+ Gy (re St) [bt 
L, 
-hd6 =o 


2) 


It is possible also to eliminate all variables except n, J and their 


derivatives from equations (2) and (4). Solving equation (4) for d@ and 


2 
d@ and substituting these into equation (2) yields: 


C Ge a hi<ew = fs 
ae Aw + (Cm, + Cm ig een Aad x 


* Cs (/- Cog JS + ACms od =: () 














Ce C Ye (9) 
Solving equation (6) for dol » and ae in terms of nand g yields: 
> ba @\. ae Cm 5 
as Gre an a . 2 (10) 
do =~ An 4 Cm 5 Ae 
Cu CLUE (11) 


Substituting equations (6), (10), and (11) into equation (9) yields: 
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z Ge Cy GC 
@ z =e oe 
dn + 8 (a me Ene (Gm," «Cn Jol dd 


Cmf a Ge m dé ia 7 


Cm oy Cras CO) 
= —_s a8 od _ 
( ir oe ) =O 


(12) 


Equation (12) if integrated once may be expressed as: 


Ci, (Ses. 29) fr = Gms (e+ “ae,) fe =v (13) 


Lx 


where 
w= Com Jy Cndg)- Cmsh yf ¢ aS (Cn +t Gn Sa oa 
( “ye C, & * Se (on 2s) a n 


Throughout the remainder of this report the following quantities 


are defined as: 











as G( =: aes | (15) 
a Cmyg + Cmae = 
Cc = Cog (/ 4 <5) (17) 
Lt = Gis = & < Ze) (13) 


Substituting equation (15) through (19) into equations (7), (3) and (13) 


yields: 
An +Bd6+CS=hbhdA (20) 
A In + BO + C/5 = hoe *, 


L/r + MA [JS = 7 (22) 


Equation (22) may be further rearranged as: 


J 
a ae | (23) 


Solution of Equations: 

In order to solve the equations of motion for values of the deriva- 
tives, one quantity must be known. In this case the known quantity was the 
eae moment of inertia about the Y axis. To determine the moment of 


inertia, the airplane was oscillated about the jack points as shown on Fig. l. 
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If it is assumed that the airplane as suspended represents a second order 
system with zero damping, it is possible to determine the moment of inertia 
about the jack points, from the formula shown on Fig. 1. The moment of 
inertia was then transferred to the center of gravity. The position of 
the center of gravity was determined by using scales under each jack point 
and under the tail for two different inclinations of the airplane. With 
the manent of inertia about the center of gravity known it was possible to 
determine fi about the center of gravity. By keeping an accurate record 
of weight added to and removed from the airplane A, was determined for each 
run. The accurate cetermination of this quantity for each run is very 
important as small errors in A cause large errors in the solution for the 
constants in the equations of motion, equations (20), (21), or (23). 
Examination of the equations of motion shows it would be most 
desirable to measure Se , « ,0/O , and SO directly. Since an angular 
accelerometer was not available, direct measurement of ae was impossible. 
Equation (6) shows the desirability for having normal acceleration data in 
order to determine Gp . For this reason, and realizing that there are 
certain inherent difficulties associated with the accurate determination of 
ad, it was decided to use normal acceleration as one of the measured 
variables. In order to make the same instrumentation system available for 
use with steady state data and so that the elevator hinge moment equation 
might be considered in future studies, it was decided to include stick 
force as one of the measured variables. The resulting instrumentation 
yielded simultaneous values of a, ny Se » a6, and Bs The actual circuits 


used are described in Appendix A. 
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Since h is a measured quantity, measurement of Bees d6, and 
ao at three different values of time would permit the solution of equation 
(20).for A, B, and C. In like manner measurement of Sn, 8, {Se and ¢cQ 
at three values of time would permit the solution of equation (21) for 
A, B, and C. Knowledge of Jn and Is » and v for two values of time 


would permit the solution of equation (22) for L and M. 


Methods of handling redundant datas 

If the data were perfect, evaluation of the variables at only 
three instants of time would be sufficient to uniquely determine the 
coefficients of equations (20) and (21). ‘With less than perfect data, 
however, some method must be used to determine the set of coefficients which 
most logically fit the data. The authors have used several methods of 
handling the redundancy arising from evaluation of the variables at many 
values of time. The first was a simple application of least squares assuming 
the error to be all in one variables the second was a method of averages: 
the third was a graphical plotting method. This latter method is the most 
desirable since it gives a physical picture of the data scatter and the 
fairing process. The lack of physical feeling for the fairing process 
when least squares or averages are used to reduce the data is of great 
importance. | 

When the method of least squares was used the normal equations of 


L 
equation (20) assuming all errors in d © are: 


AEN + BPERdO + CEnS = hEnd@ (24) 
ASind9+ BE&WAY + CE5d0- 4 Edbd9 i) 
(26) 


ASaS + BESd9 tCE SS thE S dB 





EQe 


In applying the method of least squares it was necessary to carry as many 
places as vossible in combining the variables and solving the simultaneous 
equations. No justification is offered for the apparent inconsistency 
of carrying eight to ten significant figures in the summations of the 
various combinations of variables when the variables themselves cannot be 
taken from data sheets with more than three significant figures. It is 
merely stated as a finding that as many places as possible must be used 
in order to get reasonable answers using least squares procedures. It 
follows that least squares procedures should not be attempted without the 
use of an eight to ten place calculating machine. It is also to be noted 
as a corollary of the above that least squares procecures are very unfor- 
giving of small errors in data transcription, summations, etc. 

The method of averages consists of simply using the average values 
of n, c@, 0, and a“e taken from random groupings of the individual data. 
The best results were obtained by simply grouping the data into approximately 
even thirds. The first of the three simultaneous equations then would be 
mace up of the summations of the variables over the first third of the time 
traces considered appropriate. This method consisted of dividing the data 
into three such groups. The simplicity and speed are such great advantages 
over even the simplest application of least squares that its apolicability 


shoulde be at least evaluated in any extensive data reduction problems. 


Practical Considerations 
Data for eight runs were reduced using any or several of the 
methods cescribed above. The results were consistent for different methods 


applied to the same run but the only answer which seemed appropriate was 





the value for B = (Cmsg + Cages). 


magnitude and inconsistent from run to run. 
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The values of A anc C were of improper 


With this in mind the equations 


for the unknowns A, B, C in determinant form were examined. 








F a'g =d6 =6 
4\fde fde 4d 
eg zd 
= 27 
4 Zr 2=d4 2d ( ) 
4 > 0s 
B/S) 2 /é od 
2Nn Z2d‘0 Zo 
hs ZN 2d 6 Z0 
oT) 2d'e 2d 
B = (28) 
2” Z dea ro 
Pde 7de Wd 
2 Yde 4 
In Zde 2d 
A | GA “eee Fae 
21 2c68 Pde 23) 
C = - ( 
2n ¥de 2a 
2” 206 23 
270 de Ye 
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Examination of equations (27), (23), and (29) shows that equation 
(28) is the only one where an error in the pitch rate gyro calibration 
would not affect the results. 

An examination of the pitch rate gyro revealed that considerable 
friction had developed in the bearings and it was very doubtful if the 
gyro wheel was getting up to governing speed. 

If it is assumed that the gyro maintains a constant speed during 
the 1.5 seconds during which data is being taken for any particular run, 
the reduced values of B will be correct. Values of A anc C will of course 
be directly dependent on the gyro speed for that particular rim. 

Because of time limitations it was impossible to repair the gyro, 
recalibrate, and take new data. This necessitated re-examination of the 
equations of motion in order to eliminate pitch rate as an unknown. The 
method is very similar to that proposed in Ref. 2 and called the derivative 
method. 

An examination of equation (14) shows that if (Cx, +C,., ) and 
C,.. are known, an estimation of Coms will define the quantity v . 
Equation (23) is then recognized - that of a straight line where the 
variables are lap and lof and the intercepts with the coordinate axes 
determine values of L and M. At each instant of time it is only necessary 
to compute fx ,» {(S , andzy s make a plot of equation (13) on linear 
coordinates; and fair the best straight line through the points to determine 
the intercepts. These graphical plots were extremely enlightening. They 
indicated that L can be determinec to much greater accuracy than can M. 


A blind application of least squares or averages could never indicate this 





pe 


nh 
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condition. Knowledge of this condition would be valuable for deciding 
which variable to minimize the error for application of least squares. 

It should also be noted that the answers are very insensitive to 
the estimated values used since the known value of (Cmaat C im ao) makes up 
the major part of v. 

In the process of data reduction it was attempted, with poor 
results, to use a cata instead of n for solving the moment equation. The 
fifteen hundred series runs yielded very poor results. This was accounted 
for by the fact that an accurate in-flight calibration of the angle of 
attack indicator to account for vane position error effects was never 
attained for these runs. The nature of the difficulty is described in 
Apvendix A. 

The a data was used in equation (6) to obtain values of C, + The 
use of a data necessitates taking into account the effect of up waa and 
the effect of pitching velocity. The effect of up wash is determined from 
steady flight calibration of the angle of attack indicator. The effect 
of pitching velocity is determined by measuring the distance of the vane 
ahead of the center of gravity and computing the velocity component at the 
vane due to pitch rate. Another effect that it was not possible to account 
for quantitatively was the effect of acceleration in pitch. Due to flex- 
ibility of the boom angular acceleration of the airovlane caused bending of 


the boom in such a manner that individual values of the Cr determined at 
a 


each instant of time yielded a plot as shown below. Better stiffening of 
the boom shoulc eliminate this source of error. 


As the dynamic response of the angle of attack vane was not deter- 
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mined. any natural dynamic response of the vane consisted of an error in 
the a data of undetermined magnitude. 
For the above reasons a data was used only to determine Cy from 


a 
the sixteen hundred runs. 


| we Tyme ot PNQAKI IMB 


anjulay acceleration. 






Cue 


The solution of the equations of motion by any of “We techniques 
used does not allow solution for unique values of the derivatives Coa, 
Cy C'ms » and Om 46 - In order that the results may contain 
approximate values for the above derivatives it was assumed that Cn, ey Cavs, 
and = = .5. It is important to realize that without the above assumption 
dynamic flight test procedures as described herein do not yield unique 


answers for the individual derivatives. 
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PROCEDURE 


The procedure followed in this report may be considered in the 
following parts: 1. Instrumentation, 2. Calibration, 53. Flight test, 


4. Data reduction. 


1. Instrumentation 

The design, construction, and installation of the instruments 
used in this report took approximately four months. The circuits actually 
used as well as the design criteria may be found in Appendix I. The 
reasons for the choice of variables to be measured may be found in the 
Theory and Analysis section of this report. 
ee Calibration 

The calibration curves for each circuit used may also be found in 
Appendix A. A brief explanation of inherent difficulties in measuring 
certain variables may be found in the Theory and Analysis section of this 
report. Calibration was, of course, carried out continuously as the flight 
test progressed to insure that the calibration curves remained applicable. 
The one exception to the above statement was the calibration of the pitch 
rate gyro. It was calibrated in the final stages of instrument installation 
at the Naval Air Test Center, Patuxent River, Maryland, and no further 
calibrations were performed. The difficulties encountered in the early 
stages of data reduction emphasized the desirability of continuous calibra- 
tion checking. 
3. Flight test 


The flight procedure for transient dynamic testing is extremely 
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simple and represents the greatest advantage of dynamic testing over 
steady state testing methods. Smooth air is scone in order that a 
good steady state flight condition may be attained. After attainment of 

a gooc steady state condition, the steady state is disturbed in some pre- 
determined manner. Pulse type elevator forcing functions were applied by 
simply moving the yoke a predetermined amount, holding it in the new posi- 
tion for approximately one second, and then returning it to the original 
position. The resulting pulses were poor approximations to an analytical 
pulse. However, in an equations of motion technique this was not imvortant. 
The best method founc for applying a step elevator input consisted of 
holding the stick against the force caused by a preset value of trim tab 
and then simply letting go of the stick. The elevator is so heavily 
damped that the resulting elevator motion was very nearly a step. 

All data was taken at 19.5 inches manifold pressure and about 1850 
engine revolutions per minute. Low engine speeds were required to reduce 
angle of attack boom resonance. Variations in speed were obtained by 
varying the rate of descent. Data was taken for two center of gravity 
vositions and two speeds using both pulse and step elevator forcing functions. 
4. Data reduction 

The various methods of data reduction are described in the Theory 


anc Analysis section of this report. 
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SAMPLE CALCULATIONS 


A samole of each type of calculation is carried out in this section 
to further clarify the procedure usec. The calculations are carried out 
on run #1640. The values of the variables taken from run #1640 may be 
found in Table II. Table I contains the various constants required in the 
calculations ance the calibration factors taken from the calibration curves. 

Equation (21) is 

A/r + 8e+ C/b = hee 

The values for /n > & » [Sand dmay be obtained for each voint indicated 


in Table II since: 


fo IY Ee e [ndt = o nodt (pe, seconds) 


where GR = calibration factor for n in gs/in- 


[So% sf [5 Ht (inch secomene 
[40d % = 8 = Go [bt Conch seconds) 


dao = Sa 
AY, o- = C Cys (inches) 


With values of the variables at each point three simultaneous 
equations may be obtained as follows where the summations are over the 


various points of time: 





ae 
A $s dt 4 : z ; 
pe Se eC Pa ay 2 ve 
Az frdte + Bo + CElSd% = hz dé 


/ / /é 
AE fod %e + 8 $e + CE/ddI% = hEde 


Actual calculations yielced the following simultaneous equations: 


- 2000091 A + .017371 B - .020173 C = .010376 
- .198996 A + .110355 B =- .074817 C = .020030 
- .6013869 A + .204940 B - .125027 C = .014376 


Solution of these equations simultaneously for B gives B = = .141841. 
It is clear that equation (7) as it stands could be solved using 
values of Jd, n, dO, and ane These variables in inches and seconcs units 


are also includee in Table II. 


If equation (6) is rearranged it is possible to solve for Cc, - 


a 
Cae -Qn+ Cag 
Ly = A~A/¢ 
xX 
Values of n, a, and c may be obtained from Table II using the calibration 


factors of Table I. The value of & so obtained should be further cor- 


rected for the error introduced in the vane indicator by pitching velocity. 


a 
te Le 
_. AG 

1 Gt 


—  Kelative Wind 


L 
Auk = ae 
Y 





ce 
where VY is true velocity in ft/sec. 


It is necessary to estimate Cx ¢ : Cm was estimatec theoretically to 
be -1.5. An incorrect value of C,,5 will of course be reflected in an 


incorrect value of Cy but the percentage error in C will be muck less 
a a 


than the percentage error in CHES - Using all points of Table II the 


following equation was obtained: 





oc 





ee ao S% ae Cans es 


— L/C 
Ss x 


; 


Co. ~ C14) (-3.44-79)4 LDC ES) 
a 


——-—-- 


Ot ete, oe! es 
~ 


J 7 








For this particular case an increase inC,,¢ of 50% gives an 


increase of 12% in Ci. 
a 

The terms on the right hand side of equation (13) have been esti- 
mated as shown at the end of this section for runs #1522, 1523, 1649, and 
1652. These estimated quantities have been determnined’d using values of 
(Crt CO mgs ) determined by equations of motion methoe and as rationalized 


in Fig. &. An average value of C74 anc the theoretical value of Gea 


have also been used in these estimations. Tables III, IV, V, and VI contain 
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the calculations for v ane the variables Side ana SSy - These variabdles 


are plotted and the values of L and M determined on Figs. 3 and 4. 
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RESULTS 


The following is a tabelstion of x44C most /mpor tant resel?s. 
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RESULTS 


The values of C were consistent with theoretical considerations 
and no trends were indicated as speed anc center of gravity position were 
changed. 

The damping derivatives were found with gooe consistency at the 
same speed but an increase in damping derivatives ( Ont Cm.) was founc 
_with increase in speed. 


The values of stick fixed maneuver margin oe ze mo} 
ve 


indicated a change of the proper amount when center of ee position 

was changed at constant speed. However, a decrease in stick fixed maneuver 
margin was found as speed was increased. The stick fixed maneuver potnt 
for 90 mph was 58%, for 110 mph it was 52%. 

The value of on was found with some consistency at each speed, 
but a decrease was indicated in Cy as speed was increased. The values 
of Cm, were somewhat hich in with theoretical values. 

The only immediately apparent enigmas in the results then seemed 
to arise as a result of speed changes. This would lead one to suspect an 
error in some or all of the constants used in data reduction which were 
functions of speed. These constants are 7 , and C The values of / and 
Ce used were very Peery examined and no errors could be found. 

During the in-flight calibration of the anyle of attack incicator 
a condition was found which might explain some of these changes with speed. 
A sharp increase in up wash at the wing tio was note? as speed was decreased 


below 85 mph indicated. The test aircraft was equipped with a spoiler near 


the wing root. The vossibility existed, then, that non-linearities in flow 
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ane lift were occurring during pull ups froi the 90 mph steady state. 

If non-linearities in the flow did exist, they would be reflected 
in low values of normal acceleration. Further the actual error in the 
normal acceleration would become greater with time. The recorded values 
of normal acceleration woule be less by a progressively larger amount than 
that which would have been recorded sans non-linearities. The recorded 
values of fn would have been less cueaet by as large a percentage as 
normal acceleration itself. The error in An was not as apparent as that 
for n and fr . 

The equation from which G and stick fixed maneuver margin were 
obtained was: 

yt Jy | 

Z Js Oh ee 
The quantity v as determined in Tables III, IV, V. and VI was predominantly 
made up of n, particularly after the first few tenths of a second. The 
non-linearities discussed then, if present, would have tended to make pak 
and ed too large, the error increasing with time. Since J was a 
larger quantity than Js these errors would tend to decrease the slope 
of the lines faired into the cata in Figs. 3 and 4. This could have piven 
the increased apparent stability found at lower sveeds. Also since the 
voints at later times were used predominantly in fairing Figs. 3 ane 4, 
the increased value of S at 90 mph could also be caused by these non- 
linearities. 

An attempt to evaluate the effects of non-linearities on the values 


of the damping derivatives as determined by equations of motion method lead 
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to no eefinite conclusions. It is possibie that some of the apparent 
increase in damping derivatives with speed may have been caused by low 
values of normal acceleration at the lower sneecd. The increase in damping 
derivatives with speed indicated an increase in tail efficiency. The test 
procedure required as low an RPM as possible to reduce the angle of attack 
boon resonance. Power used was 19.5 inches,full low RPM at both speeds. 
This requiree an increase in rate of descent with increase in speed. This 
does not appear consistent with the increase in tail efficiency unless the 
tail was ooerating in a particularly unfortunate position with respect to 


wing wake at the lower speed. 
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CONCLUSIONS AND RECOMMENDATIONS 


Dynamic flight testing is in general a perfectly feasible method 
of accurately determining the stability derivatives of an airplane. 
Certain basic requirements for obtaining good results are as follows: 

1. Inclusion of pitch rate as one of the measured variables. 

2. Weight, moment of inertia, and airspeed must be known very 
accurately. 

3. Instrunentation must be very accurate without unknown instru- 
ment dynamic characteristics. 

Dynamic testing has the following advantages and disadvantages. 
Advantages: 

l. Simplicity of flight procedure. 

ee Short duration of in Piiegneeresvinc. 

S- itlon-linear effects may be included in the equations of motion. 

4. Values of the derivatives obtained are the ones the airplane 
feels in response to transient maneuvers. They are the values for velocity 
equal a constant. 

5- Close correlation exists between frequency response determinec 
from dynamic testing and auto pilot design requirements. 
Disadvantages: 

l. Instrumentation is generally very expensive. 

&. Data reduction is lonz, tedious, and subject to error. 

3. Some of the results can not be compared to static flight test 
results which are for ov" equal a constant. 


In recommencing that the work be continued the most important 
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improvements in the methocs used are listed below. 

l. Further efforts should be made to ottain accurate pitch rate 
data. 

2. A vendulum or turntable should be constructed for calibration 
of the rate gyro and accelerometer. 

3. An angular accelerometer would vernit the measurement of 3“e 
directly. 

4. Angle of attach data would be improved by boom stiffening, 
anc dynamic calibration of the vane. 

o« A frequency analysis of the airplane response might yield 


interesting results. 
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APPENDIX A 


Description of Test Equipment 


The test airplane was a standard Ryan Navion instrumented to record 
tine histories of elevator angle, angle of attack, stick force, pitch rate 
and normal acceleration. These variables were recorded with a Consolidated 
Sngineering Corvoration type 5-116-P3-14 recording oscillograph. A wiring 
diagram for the instrumentation system used is shown in Fig. 5. Schematic 
diagrans showing the operation of the individual circuits appear as Figs. 

6 and 7. 

Certain basic considerations were common to the design anc operation 
of all of the circuits. These were: 

(a) To make the equivalent resistance of each circuit be 350 ohns 
as seen from the oscillograph leads, in order to meet the manufacturer's 
recommendation for 624 of critical damping for optimum galvanometer response. 

(b) To overate the sensitivity adjustment of each circuit at 
about 606 attenuation. 

(c) To use as small balancing resistors as possible in order to 
maintain linearity of the output, and at the same time keep the resistors 
large enough to limit the current drawn from the battery. 

(a) To obtain maximum galvanoneter deflection for the range of 
variables measured at each flight condition. The galvanometer deflection 
was one inch per 12.4 & amps of current. 

Power for the instruments was obtained from a 45 volt dry cell battery 
with a center tap. Voltage from this battery to the bus lines was regulated 


to 20 volts for all tests using a variable resistor and voltmeter. 
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Elevator Angle and Anele of Attack Circuits. The b. and «circuits 


ese SS —_= eth mes ——— = 








are shown schematically in Fig. 7. The potentiometers used to sense the 
angular motions were 12,500 ohm, Micro-Max precision potentiometers with 
outouts guaranteed to be linear within one per cent. The motion of the « 
vane and the elevator was transmitted to the potentiometer shafts in a one 
to one ratio. The 2500 om balancing potentiometer is used to adjust the 
circuit output to zero at various elevator or a-vane positions. The addi- 
‘tional 5009 ohm resistors in series with the balancing potentiometer were 
required to limit the current drawn from the dry cell. The 1500 ohm load 
resistor was required to adjust the current to the galvanometer to the 
desired magnitude. Sensitivity adjustment was achieved through a 500 olm 
impedance pad. The 1000 olm variable shwmt resistor was placed in parallel 
with the circuit to make the galvanometer external Camping resistance be 
550 ohms as requirec. 

The calibration of the elevator angle circuit is shown in Figs. 8 
and 9. These calibrations were made using a propeller protractor mounted 
on the elevator to measure A b.3 . Fine adjustments of the elevator 
position were made by moving the elevator with a hydraulic jack under the 
trailing edge. 

Calibration of the angle of attack circuit was made in flight using 
the proveller protractor to measure AXs, Calibrations were also made 
on the ground using a specially made fitting on the boom to measure the 
angle between the vane position and the boom centerline. These calibrations 
are shown in Figs. 10 and ll, and illustrate the difference in calibration 


under flight conditions caused by airloads on the vane and boom. Due to 





-3l- 


an oscillograph breakdown it was impossible to get an in-flight calibration 
for the 1500 series runs before the sensitivity and zero balance settings 
were changed. Therefore the in-flight calibration for the 1500 series runs 
was mace using the difference between ground and flight calibrations for 
the 1600 series runs and extrapolating for a correction to be added to the 
1500 series ground calibration. 

Stick Force, Pitch Rate and Normal Acceleration Circuits. The oper- 
ation of the Ps 8, and n circuits are similar and are shown schematically 
in Fig. 6. The strain gauge bridge circuits used to measure these variables 
will be described individually later. The 5990 or 50,900 ohm balancing 
potentiometer in each circuit gave a zero output adjustment to the circuit 
over a wide range of values of Pas 6, and n. The 5000 or 50,000 ohm re- 
sistors in series with the balancing potentiometers were necessary to vrovide 
a current limit to protect the galvanometer in case the circuit should 
inadvertently be turnec on with the sensitivity adjustment full open and 
the balancinz potentioneter turned against the stops in either direction. 

The sensitivity adjustment and damping resistors are the same as described 
for the a and §. circuits. 

The accelerometer used in the n circuit was a model F-1.5-880 accel- 
erometer manufactured by the Statham Laboratories, Incorosorated. The manu- 
facturer's specifications required an 18 volt input to the accelerometer 
which was obtained with a 100 ohm variable resistor in series with the 
circuit leads to the power supply. The accelerometer had a range of +1.5 g's. 
The circuit was calibrated between zero and one g and assumed to be linear 


to 1.5 g. The calibration curve is shown in Fig. 12. The noise level in 
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the n circuit was extremely sensitive to airplane vibrations, particularly 
to engine rpm. In order to minimize this the T-section filter shown in 
Fic. 6 was usec to filter out vibrations above 12 cycles per second. 

The stick force bridge was made up of SR-4 strain gauges cemented 
to the spokes of a specially constructed force wheel which replaced the 
standard pilots wheel in the airplane. Calibration of the stick force 
circuit is shown in Fig. ls. 

The pitch rate circuit bridge was made up of SR-4 strain gauges 
mounted on a leaf spring which restrained the precession of a gyroscope 
when the gyroscope was subjected to an angular velocity. The SR-4 strain 
gauge bridge requires a 15 volt input which was obtained by using a 100 ohm 
variable resistor in series with the circuit leads to the power supply. 

The pitch rate gyro was calibrated on a calibrating pendulum at the Naval 
Air Test Center, Patuxent River, Md., an? subsequent recalibrations were 
not feasible. This calibration established the linearity of the circuit, 
and is shown in Fig. 14. This linearity extended through six points not 
appearing within the range of Fig. 14. A 250,000 ohm +9.1% and a 59,000 om 
+).1% precision resistor was placed in parallel with one side of the bricge 
as shown in Figs. 5 and 6. At the time of the original pendulum calibra- 
tion the bridge unbalance caused by these calibration resistors was noted 
to give a circuit output equivalent to particular values of pitch rate. 
Then with the instrument installed in the airplane a calibration could be 
made for each test flight. This was done by momentarily closing a calibra- 
tion resistor switch with the aircraft in steady flight. A galvanometer 


deflection was thus obtained for a known value of pitch rate. This point 
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along with the point of zero deflection for zero pitch rate was then used 
to establish a calibration curve for the flight. 

The electric motor in the ey roscope was powered from a e4 volt 
storage battery. It was discovered during the data recuction that the 
gyro motor had hae a progressive bearing failure, and that 6 values were 
unreliable. This failure could have been ciscovered had a calibrating 


pendulun been available for calibration after each test flight. 
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APPENDIX A 


The following is a tabulatiam of the results from 
Princeton Reports No. 231 and 232. Report No. 232 consists 
of the determination of longitudinal stability parameters, 
for a standard Ryan Navion, using theoretical calculations 
and steady state flight test: techniques. Report No. 231 
consists of determining as many as possible of these same 
parametes for the same airplane using dynamic flight test- 
ing techniques. 
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